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Abstract 

The structural properties of polymers adsorbed onto a surface have been 
widely investigated using self-consistent mean-field theories. Recently, an- 
alytical mean-field theories have been applied to study polymer adsorption 
on curved surfaces but all in the context of the ground state dominance ap- 
proximation in which the polymer chain length (TV) is essentially infinite. 
Using an expression for the free energy by Semenov, we determine leading 
order (in 1/iV) corrections due to the finiteness of the polymer chain length 
on surface tension, spontaneous curvature, and rigidity constants. 
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1 Introduction 



The influence of polymer adsorption onto surfaces has been a topic that has at- 
tracted much attention for several decades. This is not only because of its great 
practical interest, but also because it serves as a testing ground of theoretical mod- 
els for confined polymer systems PUI21IIII- Recently, investigations have been dedi- 
cated to study the interplay between colloidal particles and polymers [HIqTIqTHIHIIU] ■ 
In these systems the curvature of the surface of the colloidal particle becomes a 
factor, and work has focused on the situation where either the colloidal particle 
is much larger or much smaller than the polymer coil 0J|3]. When the radius of 
curvature is large, it is reasonable to expand the free energy for polymer adsorp- 
tion in curvature. Helfrich supplied the general form of the surface free energy 
expanded to second order in the curvature [TU] : 
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where R\ and R2 are the local radii of curvature. Apart from the surface tension, 
a, one can identify three parameters in the Helfrich free energy that describe the 
physical properties of the curved interface: the spontaneous curvature 1/Rq, the 
bending rigidity constant k, and the Gaussian rigidity k. 

A great deal of attention has turned to the determination of the value of these 
curvature parameters in complex systems, both theoretically and experimentally. 
For surfaces interacting with polymers, the curvature parameters were calculated 
by a number of groups [c[10IElEl[Hl[I2J[iniEH[IS] . In the work by Eisenriegler and 
others jHllZllHlEl, the polymer density is assumed to be zero at the surface of the 
colloidal particle so that a depletion layer exists around each particle. For sur- 
faces with enhanced polymer adsorption, calculations ^21E1EI] were performed 
in the context of the ground state dominance approximation m which 

the polymer chain length is essentially infinite. In the present work, we use recent 
extensions to ground state dominance to determine corrections to the curvature 
parameters due to the finite length of the polymer chain. Such a calculation is of 
interest since the polymer chain length is an important parameter in experiments 
and computer simulations JB] thus providing a more stringent testing of theoret- 
ical models. Moreover, it was expected [211111, and later verified [201121112211221121, 



that in certain situations, the 'tails' of the polymer chain become important, lead- 
ing to qualitatively different behavior. An example is the repulsion between two 
planar walls when the separation is of the order of the polymer's radius of gyra- 
tion [22] • In the ground state dominance model the radius of gyration is infinite 
and the interaction between plates is attractive for all separations [221123 ■ 

The first extension of the ground state dominance model was proposed by 
Semenov, Bonet-Avalos, Johner, and Joanny [20] • They took the presence of 
tails into account by including a second order parameter related to the end seg- 
ment density. Good agreement was obtained for the loop and tail distribution 
of adsorbed polymer when the theoretical predictions are compared with lattice 
self-consistent mean-field calculations 21J. The two-order parameter model was 
later extended to also take into account the presence of free (non-adsorbed) poly- 
mer away from the surface [2UI2II1I2I] • Subsequently, Semenov showed that the 
two-order parameter model may be cast into a free energy formalism of a sin- 
gle order parameter [2S] in which the Euler-Lagrange equations are the Edwards 
equations, just as de Gennes had previously done in the context of the ground 
state dominance approximation dj. In this article, we use the expression for the 
free energy by Semenov [23 to determine the polymer length corrections to the 
curvature parameters. 

The outline of the paper is as follows: In the next section we review the de- 
scription of finite chain length corrections to polymer adsorption on a planar wall 
in the two-order parameter model. We show how from the two-order parame- 
ter model a single order parameter model, to which we refer to as the Semenov 
model, is constructed to describe leading order corrections in 1/N to ground state 
dominance. In section 3, the single order parameter model is applied to the study 
of curved surfaces. Explicit expressions for the surface tension and curvature 
parameters are given. We end with a discussion of results. 
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2 Semenov model for adsorption onto a planar 
surface 



Before addressing the properties of curved surfaces, we discuss finite chain length 
corrections to the structure and tension of polymer adsorption onto a planar sur- 
face. We then generalize the description to polymer adsorption onto curved sur- 
faces. 

The situation under consideration is that of a polymer chain adsorbed onto 
a planar wall located at z = (see Figure 1). The calculations presented here 
are performed in the context of self-consistent mean-field theory whereby a sin- 
gle polymer chain is considered with the effect of the other polymers taken into 
account by the presence of an external field proportional to the local segment den- 
sity [IE]- Furthermore, the interaction of the polymer with the wall is taken into 
account through an infinitely short-ranged (attractive) interaction potential pQ. 

Before taking the finite length of the polymer chain into account, we remind 
ourselves of the de Gennes-Lifshitz description P~|E| f° r polymer adsorption of 
an infinite chain, the so-called ground state dominance approximation. The free 
energy is a functional of the polymer segment density <j>{z) 0: 
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where &b is Boltzmann's constant, T is the temperature, A the surface area, and 
a is the polymer segment length. The subscripts b and w refer to the value of the 
polymer segment density in the bulk, <p b = lim <f>(z), and at the wall, <fi w = 0(0), 

2— >00 

respectively. The first term is the Lifshitz expression for the chain entropy |17j . 
the second term gives the mean-field interaction between polymer segments in a 
good solvent (v is the excluded volume parameter), and the last term accounts 
for the polymer interaction with the surface pQ (d is the extrapolation length; its 
inverse measures the interaction strength with the surface). 

The above free energy is the excess or surface free energy since it is constructed 
such that when the polymer segment density becomes equal to the bulk segment 
density, <p(z) — > 4>b, the integrand reduces to zero. The result is that the surface 
tension is equal to the above free energy per area when evaluated at its minimum, 



i.e. when the 'optimum' segment density profile, determined by the Euler-Lagrange 
equation, is inserted into it [26j. 

Rather than the segment concentration, it is convenient to consider the free 
energy as a functional of ip(z) = (p(z) l l 2 : 

"ST*- (2 ' 2) 



where ip w = ip(Q). Minimization of the above free energy gives the following Euler- 
Lagrange equation for i/j(z) with boundary condition: 

ji(>"(z) = v[^{zf-cj> b ]^{z) : 

< = (2.3) 

The Euler-Lagrange equation can be recognized as the Edwards equation for an 
infinitely long polymer chain [Ullbj. 

The expression for the free energy in eq !2.2l has been the starting point of many 
density functional treatments of polymer solutions 012] • Owing to the nature of 
the ground state dominance approximation, this treatment is, however, limited 
to infinitely long chains. Recent investigations [SHI El ESS El] have therefore 
been addressed to include finite length corrections to the ground state dominance 
model. 

In order to account for the finiteness of the polymer chain length, a second 
order parameter, <p(z), was introduced (20] connected to the description of chain 
ends. To understand better the physical origin of this term, consider a polymer 
chain with at least one segment touching the surface. Such a chain consists of a 
series of 'loops' with a 'tail' on both sides |2] (see Figure 1). The segment density 
of adsorbed chains, </>(z) is then a sum of two terms; </>i(z) = ip(z) 2 connected 
to the loop segment density and 4>t(z) = ip{z) <p(z) connected to the tail segment 
density: 

</>{z) = <j> t {z) + = ^{zf + i>(z) <p(z) . (2.4) 

This description is restricted to polymer chains that have one or more segments 
touching the surface. As a consequence, since the polymer chains are finite, 
lim <j)(z) = 0. The segment density <fi(z) therefore does not describe the poly- 
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mers in the bulk region. 
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The density of chain ends, <f> e (z), is proportional to a single ip( z ) [20: 

(j) e {z)=B^{z). (2.5) 

The value of the proportionality constant B is obtained by a normalization onto 
the number of end-segments. To show this in more detail, we denote V as the 
total adsorption of chains having one or more segments touching the wall: 

oo 

T = Jdz[^{zf + ij>(z)<p(z)]. (2.6) 



The proportionality constant B is then determined by the fact that every adsorbed 
chain has two ends so that the fraction of end-segments therefore is 2/N, with N 
the total number of segments (2Hj: 

oo oo 

JdzB^{z) = ^T = ^ Jdz[^(z) 2 + ^(z)cp(z)]. (2.7) 



The profiles of the two order parameters, ip(z) and (p(z), are determined by 
the following set of Edwards-like equations 
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— ^"{z) = v[^{zf + i>{z)tp{z) + e-<f> b }il>{z) 
o 

= v[i:(z) 2 + i:(z)^z) + e-^] V (z)-B. (2.8) 

o 

As boundary conditions at the surface, we have: 

V4 = -^, ^ = 0. (2.9) 

The latter condition follows from the interpretation of ip(z) <p(z) as the tail seg- 
ment density distribution, that is zero at the surface by definition. As boundary 
conditions far away from the surface, we have: 

lim ipiz) = , lim (p(z) = — — . (2.10) 

z^oo z^oo Vie — (pb) 

The following functional is derived as a first integral of the two differential equa- 
tions in eq 12.81 identified as the free energy: 
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(For notational brevity we suppress the ^-dependence of ip, and ip in the free 
energy.) 

The parameter e, which appears in the differential equations for ip(z) and (p(z), 
may be interpreted as the chemical potential (in units of k^T) of the attached 
polymers. Its value is determined by a balance of the chemical potential of the 
polymers at the surface with the chemical potential of bulk polymer. It is derived 
in ref. |2*m i2*T] that e is approximately given by: 



o 

Combining the above expression for e with eqEZl we can rewrite the normalization 
condition for B as: 



This model is the two-order parameter model which has been shown to be in 
agreement with numerical solutions of the self-consistent field equations [21] . One 
disadvantage of the two-order parameter model is the rather artificial distinction 
made between chains that have one or more segments touching the wall and free 
polymer chains, leading to the presence of e^O. Furthermore, the above formalism 
to determine the segment density profiles and free energy does not correspond to 
a free energy functional formalism in which the differential equations determining 
ip{z) and <f(z) are the Euler-Lagrange equations to the free energy in eq 12.111 The 
reason is that a minimization of F would not necessarily lead to the boundary 
condition ip(0) = 0, and, more importantly, B and e are not constants in the 
minimization but are themselves functionals oi ip{z) and ip(z) (eqs 12.71 and 12.121) . 

To circumvent these difficulties, we now turn to the derivation of what we term 
the Semenov model j2H|. In this model, no explicit distinction between adsorbed 
chains and free chains is made leading to a single order parameter describing the 
total segment density, <p(z). As a first step we consider again the free energy given 
by ea l2~TT1 
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Since adsorbed and free chains are now treated equally, we have taken e = 0. 
Furthermore, we have set = 0, for the time being, and deal with the inclusion 
of bulk polymer below. 

Next, since the combination ip(z) (p(z) is related to the tail segment density, 
whose contribution to the total segment density profile is of order 1/N, we can 
conclude that ip(z) describes effects that are of order 1/N, i.e. ip(z) = 0(1/N). 
With (f>(z) = ip(z) 2 + i>{z) the total segment density, we can therefore write 
in an expansion in 1/N: 
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Inserting this into the free energy in eq 12.141 we have 
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with corrections of 0(1/N 2 ). Next, we insert the explicit formula for B given in 
eq 12.131 with e = 0, into the free energy. We are then left with the following free 
energy functional: 



Ak B T 



dz 



a 2 W 
24 6 



G(<f>) 



1^ 
d~6 



where 
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(2.18) 



To impose a certain value of the bulk polymer density, we introduce the bulk 
chemical potential in G((f>): 



G(<P) = -ct> 



2_ ,1/2 ,l/2_ 



(2.19) 



2 r N Vb 

The value of the bulk chemical potential, /x, is chosen such that the bulk density, 
as given by the minimum of G(<ft), is equal to a certain 0^: 



G'(0 b ) = ^ - ^ - /i = 



(2.20) 



As expected, the chemical potential is lowered by reducing the chain length. As a 
final step, we subtract from G((f>) the asymptotic bulk free energy G(<fib) in order 
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for the above free energy to be the surface free energy. Our final expression for 
G(<f>) then reads: 

G(0) = |(0-0,) 2 + ^(0 1/2 -0; /2 ) 2 . (2.21) 

One notices that the above form oiG(4>) reduces to that in ground state dominance 
(See eq 12. 1|) in the limit N^oo, as it should. 

Even though we started from an expression for the free energy in eq 12.111 
which we stressed is not the free energy functional, we have now constructed a 
free energy functional in eq 12.171 It is a functional of <j)(z), that is now the total 
segment density profile in which no reference to 'loops' or 'tails' is explicitly made. 
Its form was first proposed by Semenov j2H] who also determined the next order 
contribution and showed it to be of 0(1/ TV 3 / 2 ). Semenov j2H] points out that this 
approach agrees with the two-order parameter model (to 0(1/N)), even though 
this agreement has to be interpreted somewhat loosely since a direct comparison 
is difficult to make owing to the difference in the treatment of bulk polymer. 

It is convenient to rescale all densities by 4> b and all lengths by the bulk corre- 
lation length ^ b = a/^3 v (fib- 

X EE Zji h , d EE d/& , (j){x) EE (j) b f (x) 2 . (2.22) 

Furthermore, we introduce as a small parameter the (square of) the bulk cor- 
relation length, £ 6 , divided by Rq ee y / Na 2 /Q, the polymer's radius of gyration: 



In order for e to be small, the bulk polymer density v<f> b '^>2/N, i.e. the expansion 
is expected to fail for very dilute polymer solutions. For concentrated polymer 
solutions, where the bulk polymer density is much higher than the overlap con- 
centration, v<p b ^> 1/ VN, this type of mean-field theory is also expected to break 
down. Our model is therefore most relevant for semi- dilute polymer solutions. 
The free energy in eq 12.171 is written as: 

oo 

FIM = ^4H= [dxU) 2 + 9(f )}-lf 2 o, w , (2.24) 
vt, b (pi Ak B T J d 
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where g(fo) is given by 

s(/ ) = (/ 2 -l) 2 + e(/ -l) 2 . (2.25) 
The Euler-Lagrange equation to the free energy in eq 12.241 reads 

f'^) = \g'{k) =► ax? = g{h), (2.26) 

which gives 

fo(x) = -VgUo) = -(/o - 1) (/o 2 + 2/o + 1 + ^) 1/2 ■ (2.27) 

As boundary condition to the above first order differential equation, we have that 
fow = ~fo,w/d- The initial value f 0)W = /q(0) is thus obtained by solving the 
following algebraic equation 

-~fo, w = (fo, w ~ 1) (fo, w + 2f , w + 1 + e) 1 ' 2 . (2.28) 
a 

The profile obtained as an explicit solution to the differential equation reads: 

^ ^ ~^ 16 exp(/3(x + x w )) — 16 — e exp(— (3{x + x w )) ' 
where we have defined 

, a v 2 + 2/ 0jto + e + a{3 
exp([3x w ) = — , 

4 [J0,w — J- J 

a = (/o% + 2/o,«, + l + £) 1/2 , 

/? = (£ + 4) 1/2 . (2.30) 

This analytical result for the order parameter profile was derived without making 
any further approximations. However, it should be kept in mind that the free 
energy in eq 12.171 captures only the leading order correction (in 1/N) to ground 
state dominance. The result is that only the two leading terms in an expansion 
in £ are physically relevant: fo(%) = fo,o(%) + £ fo,i( x ) + • • •> with 

/o,oG«0 = I w 1 , x , (2-31) 
tann(x + xq) 

(5 e 4x ° - 8 e 2x ° + 4 - e~ 4x °) (4x + 4 e" 2a '- 2a;o - e~ Ax - Ax ° - 3) 
~ ~ 32 (e 4 *° + 1) sinh 2 (x + x ) 32 sinh 2 (x + x ) ' 
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where we have defined Xq = (1/2) arcsinh(2<i). These expressions can be used to 
analyze the limiting behavior of the segment density profile for strong and weak 
adsorption. For weak adsorption, we have 



[Z) 1 + 4 



e -2z/& _ £ (1 + H£) e -a*/& + 
8 £& 



(2.32) 
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In the case of strong adsorption, one finds the characteristic l/z 2 -behavior pQ of 
the segment density profile close to the solid surface 



<Pb z 2 



b 



(2.33) 



This result is unaffected by the finite chain length corrections which implies that 
for strong adsorption, the segment density profile in the immediate vicinity of the 
wall is unaffected by the length of the polymer chain. 

The surface tension is given by the surface free energy with the order parameter 
profile, as determined by the Euler-Lagrange equation, inserted into it: 

2 a 



a = 



2 Jdx[f {xf]-lfl w . (2.34) 



The surface tension can be conveniently rewritten as an integral over / 

fo,w 

~a = 2 JdfoV^(foj-ifl w . (2.35) 
i 

The value of the surface tension is thus calculated either by inserting the density 
profile in eq 12.291 back into eq 12.341 or by direct evaluation of the density integral 
formula above. In either case, the result is: 

~ a /„,„ (fo, w - 1) + ~ (2 - e) ((3 - a) - 2s In ( 1+ /+/ Q ) " ' (2 ' 36) 

This analytical result for the surface tension was derived without making any 
further approximations. However, it should again be kept in mind that only the 
two leading terms in an expansion in e are physically relevant: 

a = a + eax + ... . (2.37) 

Expanding eq 12.361 in e gives as explicit expressions: 

4 4 cosh(xo) 



a 



3 3 tanh(a; ) 3 sinh 3 (x ) ' 



di = =s6o- 1+21n < 1 - e "* , >- (2 ' 38) 
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The result for ctq was first derived in ref. in Figure 2, the surface tension is 
plotted as a function of chain length at fixed d— 1 and for various values of vfo. 
The solid line is the full result in eq 12.361 with the dashed line the two leading 
terms feqs 12.371 and |2~3%|) . Even though the relevant variable is the combination 
e = 2/(Nv (fit,), we have chosen to plot the surface tension as a function of N for 
various vfo. The condition that e is small implies that iV ^> 1/ (v (fib e) . This limit 
corresponds (for each vfo) roughly to the point where the solid and dashed lines 
deviate. The leading order correction to the surface tension is always positive 
indicating that the lowering of the surface tension due to the presence of adsorbed 
polymer is reduced when the polymer chain becomes shorter. The magnitude of 
this effect is appreciable (~ 10 %) already for large polymer chains (N ~ 1000) 
when the bulk density is low (t>0& — 0.001). 

It is instructive to consider the behavior of a± in the weak adsorption limit 
(<*»&): 

ffl = 16rf2 + 48^ + C>( rf4 ) ' (2 ' 39) 
This result for b\ can be combined with the limiting behavior for <7 . To show this 
in more detail, we subtract from the surface tension the constant contribution to o 
that remains even when the density profile is equal to the bulk density everywhere, 
(fi(z) = (fib. One can then write 

, m a 2 (fit, a 4 kpT / e \ , s 

A ^ + , Br _^ = ___!_(! __ + ...). (2 . 40) 

This expression more clearly shows that the lowering of the surface tension due 
to the presence of adsorbed polymer is reduced when the polymer chain becomes 
shorter. One may show j2jj that the next term in the expansion of the surface 
tension in eq 12.401 goes as ~ — 0.1540 e 3 ^ 2 . 

For completeness we also give the expression for a\ in the case of strong poly- 
mer adsorption (cf<c£b) : 

a x = j + (2 ln(2 d) - l) d + 0{d 2 ) . (2.41) 
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3 Semenov model for adsorption onto curved 
surfaces 



We now extend our analysis to curved surfaces. In particular, we are interested in 
determining the coefficients of an expansion of the free energy to second order in 
the curvature: the spontaneous radius of curvature, Rq, and the rigidity constants 
k and k. In general geometry, the free energy functional in the Semenov model 
reads: 



knT 



dr 



a 2 |V0(r)p 



24 
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Aa?_ 
d 6 



(3.1) 



We rescale the variables as done previously, see eq 12.221 The rescaled free energy 
is then given by: 

1 
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(3.2) 



with the Euler-Lagrange equation and boundary condition: 



1 



A/(f) = -</(/), 



n ■ V/™ = -\f u 
d 



(3.3) 



To derive the curvature parameters, both order parameter profile and free en- 
ergy are expanded to second order in curvature for a spherically and cylindrically 
shaped surface [28] . For example, in the spherical geometry we have: 



f.(x) = f (x) + ^ f Sjl (x) + — /, i2 (x) + . . . 



(3.4) 



where R is the sphere's radius. In this expansion, the Euler-Lagrange equation 
and boundary condition are given to zeroth and first order by: 



d 

~J s,l,w ■ 

d 



(3.5) 
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The spontaneous radius of curvature, Rq, and the rigidity constants k and k, are 
given by [11113123 



Co 
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k =2 [dx[x 2 f Q {xf] , 
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- jdx [fo(x)f sA (x)} 





(3.6) 



The expressions for the spontaneous radius of curvature, i?o, and the rigidity 
constants associated with Gaussian curvature, k, can be rewritten as integrals 
over /, similar to the expression for the surface tension in eq 12.341 

fw fu 



Co 



2 UtVgU) 



df 



1 

7W) 



(3.7) 



One may show that the bending rigidity constant, k, is given by 

Ad 

k 



-, 2 



dfvauj 



V / ^[2 v ^)-rfV(/»)] 

fw fw f 



(3.i 



Explicit results for the value of Co, k, and k is obtained by inserting the functional 
form for g(f) feq I2.25|) into eqs 13.71 and 13.81 and evaluating the resulting integrals 
numerically. Typical results for the various curvature parameters as a function of 
chain length are shown as the solid lines in Figure 3. Again, it should be kept in 
mind that the physical relevance of these expressions for the curvature parameters 
is restricted to the first two terms in an expansion in e: 

Co = c ,o + e c 0) i + . . . , 
k = k + eki + ... , 

k = k + e~h + . . . . (3.9) 
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An explicit calculation gives [To] : 

c 00 = \ ln(l - e~ 2xo ) , 

3 sinh (xq) 3 

% = -- + TT \, S -|dilog(l-e- 2 -°), 
3 3 tann(x J 3 

(27e 2 "° + 5 - e~ 2x ° + e" 4 " ) 
fc ° " i8(e2-o_i)( e 4xo + i) ' 



and 



5c,! = ~ V 12(eto0 + i) / +dilog(l-e-^) + -ln(l-e-^), 

| (3e 2x ° - e- 2xo ) _ (5e 4x ° - 8e 2x ° + 4 - e" 4x °) . _ g _ 2 

1 12 (e 4 ^ + 1) 12 (e 4a; o + 1) n l ~ e J 

+^dilog(l- e - 2 ^)-Li 3 (e- 2 ^), 
o 

jfei = r^- ^ [456 e Wx ° - 354 e 8x ° + 1272 e 6xo - 281 e 4xo + 792 e 2x ° 

288 (e 4a: o + l) 3 L 

-143 + 168e~ 2a '° + 37e~ 4x ° - 3e" 8l °] (3.11) 

(13e 4 *"-16e 2a;o +4-e- 4a;o ) _ 2a;o . 1 ... _ 2:Eo . 

- — : r ln(l - e ^ x °) + - dilogfl - e , 

12(e 4a; o + l) V ; 2 &v ; ' 

where the polylogs Li n (z) and dilog(z) are defined in the Appendix. 

The two leading contributions eq I3.1UI with eq 3.11 are shown as the dashed 
curves in Figure 3. Again, the validity limit of the expansion in e corresponds 
(for each v<fib) roughly to the point where the solid and dashed lines deviate. 
The curves in Figure 3 show that the leading order correction to the curvature 
parameters is significant already for large polymer chains, especially when the bulk 
polymer density is low. Furthermore, when the polymer chain becomes shorter, 
the effect of polymer adsorption on the curvature parameters is reduced. However, 
this is not a general result and it is only valid when the adsorption strength is 
sufficiently weak. This is shown in Figure 4 where we have plotted the leading 
order corrections to the curvature parameters as a function of the adsorption 
strength. For large d, i.e. weak adsorption, the leading order corrections are such 
that they are opposite in sign to the ground state dominance results (co,i and k\ 
are positive and k\ is negative); for small d, i.e. strong adsorption, they enhance 
the effect. 

These results are further demonstrated by considering the limiting behavior. 
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In the weak adsorption limit (d £&) one finds: 

11 11 <n( 1 

^ " 32^ + 288^3 + ° ( ^' 

fci = — H — + ), 

128 576 rf s V 

9 1 11 1 
fci = — 4-— + (3.12) 
256 ^2 96 d 3 v rf 4 ; v 7 

The leading order corrections all scale as 1/d 2 just as the ground state dominance 
results for £0,0, fco, and fco [T3]. We can therefore write: 

2k a 4 fc B T , £ , 

- a 4 fc B T^ 6 , _ 3_£ v 

~ 288 wd 2 1 8 

fc = -— ^(1-- + ...)- ( 3 - 13 ) 
192 vd 2 V 8 ; V 7 

These expressions show more clearly that in the weak adsorption limit, the influ- 
ence on the curvature parameters caused by the presence of adsorbed polymer is 
reduced when the polymer chain becomes shorter. 
For strong adsorption (d<^£ 6 ) one has: 

Q,,, = ? ln(2d) + j + i + (2 ln(2d) - 3) d + 0(<?) , 

27 7T 2 

h = -^ + 7^ + ( ^ 2 )> ( 3 - 14 ) 

where the zeta function ((z) is defined in the Appendix. 

These results show that when the adsorption is sufficiently strong, one finds 
the somewhat surprising result that shortening the polymer chain enhances the 
effect on the curvature parameters induced by adsorbing polymer. 



4 Discussion 

In this article we have investigated the chain length dependence of the curvature 
properties of polymer adsorbed surfaces. Our calculations are done in the context 
of the Semenov model j2H| for the free energy which captures the leading correction 
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to the free energy in an expansion in 1/N. Analytic expressions are derived for 
the leading corrections to the surface tension and curvature parameters. 

For an infinite chain, it is well-established that the ground state dominance 
model predicts that the addition of polymer reduces the surface tension, induces 
a spontaneous curvature towards the polymers, gives a positive contribution to 
the Gaussian rigidity, and reduces the value of the bending rigidity. We find that 
if the adsorption strength is sufficiently low, the leading correction due to the 
finiteness of the polymer chain to the surface tension and curvature parameters is 
to reduce the effect of polymer adsorption on these coefficients. The magnitude of 
this effect may be considerable already for long polymer chains, especially when 
the bulk polymer density is low. 

For strong adsorption, we find that the leading contribution to the curvature 
parameter changes sign. This means that shortening the polymer chains now 
enhances the effect polymers have on the curvature parameters. This seems a bit 
surprising since it is expected that for very short chains, the influence of these 
chains becomes less pronounced. Still, our result is not necessarily in contradiction 
with this expectation since for small N all the terms in the expansion in 1/N are 
important. 

One would like to understand the physical reason behind this "enhancement 
effect" . Mathematically, the effect is traced back to the non-monotonic behavior 
of the leading correction to the segment density profile fo,i(x), which is directly 
linked to the leading corrections to the spontaneous curvature (c 0j i) and the Gaus- 
sian rigidity (ki). (An equivalent argument can be made for the bending rigidity 
(ki) which is determined by the behavior of 

In figure 5, we have plotted fo,i(x) for various values of the adsorption strength. 
For weak adsorption, fo,i(x) monotonically increases to zero at large distances 
indicating that short chains adsorb less. For strong adsorption, two regions may 
be identified: in the vicinity of the wall the segment density profile is dominated by 
polymer 'loops'. In this region fo,i(x) is close to zero indicating that the polymer 
density is independent of chain length. At somewhat larger distances from the 
wall the profile is dominated by the polymer 'tails' [2, 24]. Here fo,i(x) goes 
through a minimum and then monotonically increases to zero at large distances. 
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One may therefore conclude that, although the total adsorption is less when the 
chain length is reduced, the segregation between a region dominated by loops and 
a region dominated by tails is more pronounced leading to steeper gradients in 
the polymer segment density profile. This is then ultimately responsible for an 
enhancement of the polymer contribution to the curvature parameters in the case 
of strong adsorption. 

We end with a discussion of the limitations of our theoretical treatment. Cor- 
rections to the curvature parameters were calculated to leading order in 1/N. We 
showed that the leading correction to the segment density profile gives a modi- 
fied density profile that varies on the scale of the bulk correlation length, as does 
the ground state dominance profile. One expects [231 127] higher order corrections 
(0(l/N 3 / 2 )) to give modifications to the segment density profile on the scale of 
the polymer's radius of gyration. Although these higher order corrections domi- 
nate the segment density profile when z 3> £&, they remain subdominant when it 
concerns the curvature parameters since these are expressed in terms of integrals 
over the entire region of inhomogeneity. 

A subtle point concerns the dimensional dependence of our analysis. Since 
all our results are derived within mean-field theory, it is implied that the criti- 
cal exponents are those of a polymer system embedded in 4 dimensions. In our 
derivation of the curvature parameters we explicitly consider the geometry to be 
that of a 2-dimensional surface curved in 3-dimensional space, which is arguably 
inconsistent with the mean-field assumption. It seems hard to avoid such an incon- 
sistency, however, and it is a critique which applies to all mean-field calculations 
of curvature properties. 

The interaction with the solid surface is modeled by the extrapolation length 
d giving an effective interaction. Details of the structure of the wall (e.g. surface 
roughness, surface heterogeneity) and the full wall-polymer interaction potential 
are therefore not explicitly considered. It is expected that these factors influence 
the polymer segment density in the very vicinity of the wall. In the present 
calculations, the scale of the segment density profile is set by the bulk correlation 
length. One expects that the detailed structure of the wall is unimportant as long 
as the bulk correlation length is much larger than the microscopic length scale 
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associated with this structure. 

Formally, in the calculation presented the surface is taken to be that of an 
undeformable, attractive solid wall so that the main application of our results to 
experiments is for, say, polymer adsorption onto a colloidal particle. However, 
our results should also be of interest to the description of polymer adsorption on 
more flexible surfaces such as membranes [2H]> for which the essential physics is 
captured by the Helfrich free energy. 

A Mathematical Functions 

As a reference, we provide definitions of the special functions used. 
The Polylogarithm function is defined as: 

OO 

k=l 

The Dilogarithm function is defined as: 

z 

dilog(z) = Li 2 (*) = Jdtj^r • (A.2) 
1 

The Riemann zeta function is defined as: 

C(n) ee Li n (l)EE^-, (A.3) 

k=l 

C(3) ~ 1.202056903... (A.4) 
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Figure 1: Sketch of a polymer chain adsorbed to a planar, solid surface located at 
z — 0. The polymer chain consists of a series of 'loops', with both ends adsorbed 
to the substrate, and with two 'tails', having one end free. 




Figure 2: Surface tension a as a function of chain length. The adsorption strength 
d—1; the bulk polymer density ranges from v<f)f, = 0.001 to v<pb = l. The solid line 
is eq 12.361 the dashed line is eq 12.371 combined with eq 12.381 
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Figure 3: The curvature parameters as a function of chain length; (a) is the spontaneous 
curvature cq, (b) is the Gaussian rigidity k, (c) is the bending rigidity k. The adsorption 
strength d=l; the bulk polymer density ranges from v4>b = 0.001 to v(f>b = l- The solid 
lines are numerical evaluations of the integral expressions in eas l3.7l and l3.8l The dashed 
lines are the two leading contributions, ea 13, 101 combined with eq 3.11. 
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Figure 4: Leading order corrections to the curvature parameters: 60,1, k±, and hi, 
plotted as a function of the adsorption strength (strong adsorption when d <C 1). 
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Figure 5: The order parameter profile /o,i( x ) as a function of x for various values 
of the adsorption strength; solid line: d — 0.01, dashed line: d = 0.1, dot-dashed 
line: d—0.5. 
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